Detection of the 205 [Nil] Line from the Carina Nebula 
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ABSTRACT 

We report the first detection of the 205 ftm 3 Pi — > 3 P [Nil] line from a ground-based observa- 
tory using a direct detection spectrometer. The line was detected from the Carina star formation 
region using the South Pole Imaging Fabry-Perot Interferometer (SPIFI) on the Antarctic Sub- 
millimeter Telescope and Remote Observatory (AST/RO) at South Pole. The [Nil] 205 /mi line 
strength indicates a low-density (n ~ 32 cm -3 ) ionized medium, similar to the low-density ion- 
ized halo reported previously in its [OIII] 52 and 88 /mi line emission. When compared with the 
ISO [CII] observations of this region, we find that 27% of the [CII] line emission arises from this 
low-density ionized gas, but the large majority (~ 73%) of the observed [CII] line emission arises 
from the neutral interstellar medium. This result supports and underpins prior conclusions that 
most of the observed [CII] 158 /im line emission from Galactic and extragalactic sources arises 
from the warm, dense photodissociated surfaces of molecular clouds. The detection of the [Nil] 
line demonstrates the utility of Antarctic sites for THz spectroscopy. 

Subject headings: infrared: ISM: lines and bands — HII regions: individual (Carina Nebula) — spec- 
troscopy 



1. Introduction 

The far-infrared (30 /im < A < 350 /im) con- 
tains a wide variety of fine-structure lines aris- 
ing from the ground-state configurations of astro- 
physically abundant atoms and ions. Among the 
brightest of these are the 205 /mi and 122 /mi [Nil] 
and 158 /im [CII] lines. 

Since the ionization potential of nitrogen is 
14.53 eV, N + is only found in HII regions. The 
ground-state 3 P term of the N + ion is split by the 
spin-orbit interaction into the three 3 P2,i,o lev- 
els from which the 121.898 fim ( 3 P 2 -^ 3 Pi) and 
205.178 /tin ( 3 Pi -^ 3 Po) lines arise. For electron 
temperatures of 8000 K, the critical densities are 



293 and 44 cm" 3 for the 122 and 205 /im lines, re- 
spectively, so that the line ratio is an excellent den- 
sity probe for low-density ionized gas. Both of the 
[Nil] lines were strongly detected from the Milky 
Way by the COBE FIRAS experiment. After the 
158 /im [CII] line, these lines are by far the bright- 
est lines in the 100 /im to 1 mm FIRAS band- 
pass, with line luminosities approximately 0.045% 
and 0.03% of the total Milky Way far-IR luminos- 
ity for the 122 /im and 205 /im lines, respectively 
(Wright et al. 1991). ISO observations show that 
the [Nil] 122 /im line is very bright in external 
galaxies as well, especially in the disks of normal 
galaxies (e.g., M83 and M51, Kramer et al. 2005). 



Despite the strength of the lines, the [Nil] 
lines were the last of the bright far-infrared fine- 
structure lines observed from NASA's airborne ob- 
servatories in the 1979-1995 period. This is be- 
cause the 122 /xm line lies very close to a strong 
telluric water vapor absorption feature, so that tel- 
luric transmission is poor (~ 30%) even at high 
aircraft altitudes, and because the 205 /im transi- 
tion frequency was poorly known. The 205 fim line 
lies near the long wavelength cutoff for stressed 
Ga:Ge photoconductors, strongly limiting the sen- 
sitivity of spectrometers that employ these de- 
vices. However, the 122 fim [Nil] line was first 
detected from the Kuiper Airborne Observatory in 
1989 and the 205 /im line in 1991 from the Galactic 
HII region G333.6-0.2 (see Colgan et al. 1993, and 
references therein). Since then, the 122 ^.m line 
was detected and mapped from a wide variety of 
sources using the Long Wavelength Spectrometer 
(LWS) on the Infrared Space Observatory (ISO). 
The 205 fim line remains relatively unexplored, 
since the LWS has a long wavelength cutoff of 197 
fim. 

The 205 fim [Nil] line provides a vital key to the 
interpretation of the 158 /im [CII] line. The [CII] 
line is important. It typically amounts to 0.3% 
of the observed far-IR luminosity for star-forming 
galaxies like the Milky Way (cf. Stacey et al. 1985, 
1991; Crawford et al. 1985; Wright et al. 1991), 
and is the dominant coolant for much of the neu- 
tral interstellar medium (ISM) in galaxies includ- 
ing atomic clouds, the warm neutral medium, and 
photodissociation regions (PDRs) on the surfaces 
of molecular clouds exposed to stellar UV pho- 
tons. However, modeling of the neutral gas cooling 
is complicated by the fact that a significant frac- 
tion of the observed [CII] line emission may arise 
from low-density ionized gas regions. Even for the 
Milky Way, it is uncertain how much of the ob- 
served [CII] line emission arises from low-density 
ionized gas, and how much arises from the neu- 
tral ISM. Ionized gas phase contribution estimates 
range from 10 to 30% (Stacey et al. 1985) to near 
50% (Petuchowski and Bennett 1993, Hcilcs 1994), 
depending on assumed gas-phase abundance of C + 
and ionized gas density. Fortunately, the far-IR 
fine-structure lines of N + offer a method for un- 
raveling the various contributions. 

The [Nil] 205 fim line and the [CII] 158 a™ 
line (n crit ~ 46 cm" 3 , T electron = 8000 K) have 



nearly identical critical densities for excitation in 
ionized gas regions. Therefore, the line ratio in 
ionized media is essentially only a function of the 
N + /C + abundance ratio. This is insensitive to the 
hardness of the stellar radiation fields since the 
photon energies required to ionize each species to 
the next ionization states are similar: to form N ++ 
and C ++ requires 29.6 and 24.4 eV, respectively. 
Therefore, the observed [Nil] line emission reveals 
the fraction of the observed [CII] emission that 
arises from ionized gas, thereby better describing 
the cooling of both the ionized gas and the neutral 
ISM. The [Nil] observations, when combined with 
[CII] work, therefore constrain the cooling of two 
major components of the ISM. 

Here we report the first detection of the 205 
fim line from a ground-based observatory using 
a direct detection spectrometer 1 . The line was 
detected from the radio peak of the Carina II 
HII region in the Carina Nebula using the South 
Pole Imaging Fabry-Perot Interferometer (SPIFI, 
Bradford et al. 2002) on the Antarctic Sub- 
millimeter Telescope and Remote Observatory 
(AST/RO) at South Pole (Stark et al. 2001). 
We use these observations to probe the density 
structure of the nebula and to determine the frac- 
tion of the observed [CII] line emission that arises 
from the ionized medium. 

2. Observations 

The Carina molecular cloud complex is a scries 
of star-forming molecular clouds that lie at a dis- 
tance of about 2 kpc along the Carina spiral arm of 
the Galaxy. In the optical band, the Carina Neb- 
ula (NGC 3372) is quite prominent, subtending 
about a square degree on the sky. At the core of 
the region lie the bright HII regions Carina I and 
Carina II, illuminated respectively by the Trum- 
pler 14 and Trumpler 16 star clusters. These two 
clusters contain 36 O stars, including 6 of the 17 
known 03 stars in the Galaxy, and 3 WN-A stars 
(see Davidson & Humphreys 1997, and references 
therein). The most famous member of Trumpler 
16 is the enigmatic, luminous blue variable star 
r\ Carinae. Eta Carinae is a single star with a 
mass in excess of 100 M Q ; it brightened by more 



1 There is a near simultaneous detection of this line reported 
by Kawamura et al. (2006) made with a heterodyne re- 
ciever. 
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than 5 stellar magnitudes for two decades in the 
mid-19" 1 century. With a current luminosity in 
excess of 10 6 L Q (Cox et al. 1995), even at a dis- 
tance of 2.3 kpc the dust-enshrouded star is the 
brightest extrasolar-system infrared source in the 
sky. The cluster's extraordinary concentration of 
O stars is unique within the Galaxy and rivals the 
30 Doradus region of the Large Magellanic Cloud 
in total luminosity. 

We mapped the Carina I and II regions in the 
205.178 ^m [Nil] line using SPIFI on the AST/RO 
telescope at South Pole. SPIFI was deployed on 
AST/RO for the 2005 Austral winter and was op- 
erational from 3 August until 7 September. A de- 
tailed analysis presenting all of the Carina Nebula 
[Nil] data will appear in a future paper (Oberst 
ct al., in prep.). Here we present data obtained 
at the line peak, and at a position located 1.27' 
from the peak that is the nearest spatial position 
sampled by the Infrared Space Observatory (ISO) 
satellite. The line peak data are a coaddition of 
3 separate pointings obtained on 2005 August 23 
and 24, while the ISO position data are a coad- 
dition of 5 separate pointings obtained on 2005, 
August 20, 23, and 28. 

SPIFI's detective elements are a 5 x 5 array of 
silicon bolometers held at 60 mK and arranged 
on a rectangular grid with a 65" spacing between 
pixels. The array is fed by an array of Winston 
cones yielding a 54" beam (solid angle, tlbeam = 
5.4 x 10~ 8 sr). The beam size was roughly verified 
through scans of the limb of the Moon. We cm- 
ployed a velocity resolution of 61 km s^ 1 (FWHM) 
for these observations, and scanned 7.54 resolution 
elements, Nyquist sampled. The velocity standard 
was set by observing the 205.4229 /im laser line 
used as the local oscillator source for the TREND 
laser (Yngvesson et al. 2004), also deployed on 
AST/RO for this season. We obtained 100 scans 
for each of three distinct pointings for a total of 
12.3 hours integration time on the line peak posi- 
tion. The three pointings were dithered by 2 pixels 
in the array to minimize flat-fielding errors. At the 
ISO position, we obtained 447 scans in 5 distinct 
pointings totalling 18.3 hours of integration time. 

Telescope pointing was initially achieved through 
raster scans of the Sun at low elevations in the 350 
/xm window during the Austral summer, refined 
in August through observations of the Moon, and 
finally verified by observations of the CO (7— >6) 



line in G333.6-0.2. Final pointing accuracy was 
<~ 1'. The array was flat fielded by rotating our 
chopper wheel at the entrance window of SPIFI to 
alternately sample between the chopper blade and 
the (very emissive) sky at 205 fim. Line-of-sight 
transmission was obtained by performing skydips 
of the telescope at 205 [im on 5 September and 
linking the skydips to the NRAO 350 /im tipping 
radiometer using the scaling laws of Chamberlin 
et al. (2003). We estimate the 205 /im line-of- 
sight transmission towards the Carina Nebula was 
between 4.3% and 5.3% at the time of the obser- 
vations. 

The data were obtained by using a 3-position 
chop of the AST/RO telescope tertiary mirror 
with amplitude of 30', thereby sampling sky emis- 
sion at ± 30' in azimuth with respect to the source 
position. The data were calibrated by obtaining 
the system gain (mV signal/Kelvin) on hot and 
cold loads placed in the f-cone of the receiver, then 
correcting this gain by the known efficiency of the 
telescope at 200 /im {rj te i = 51%, Stark et al. 1997) 
and the measured transparency of the sky at the 
time of the observations. 

As calibrated by our chopper wheel, over the 
course of an hour integration, the best pixels 
in SPIFI achieve an instrument noise equivalent 
power (NEP) (referred to the front of the dewar) of 
2.5 x 10~ 15 W Hz" 1 / 2 . This is about a factor of 10 
better than the best NEPs obtained with stressed 
Ge:Ga photoconductors at these wavelengths un- 
der similar backgrounds (cf. Colgan et al. 1993). 
The sensitivity ratio is consistent with the much 
better detective quantum efficiency of bolometers 
(i]boi ~ 50%) compared with stressed Ge:Ga pho- 
toconductors at these wavelengths (rjGe-.Ga ~ 3%, 
cf. Stacey et al. 1992) and the presence of 
generation-recombination noise in the photocon- 
ductors (cf. Rieke 2003). SPIFI's sensitivity cor- 
responds to a double-sideband heterodyne receiver 
temperature of ~ 150 K, or a factor of ~ 7 better 
than the best reported values at 1.4 THz (Yngves- 
son et al. 2004). 

3. Results 

We have three independent pointings during 
which a pixel was pointed at the radio contin- 
uum peak of the Carina II nebula (1 = 287. °57, 
b = -0.°59; RA =10' l 45 m 01. s 0, Dec = -59°38'13" 
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(J2000), Retalack, 1983). The spectra from these 
three pixels were individually corrected for tel- 
luric transmission and aperture efficiency, and cal- 
ibrated as described above to yield units of Tm b ■ 
The [Nil] line is clearly detected with a statistical 
significance of 20 standard deviations (Figure 1). 
The peak antenna temperature is 0.42 K, and the 
integrated line flux is 1.03 ± 0.05 x 10" 14 W m" 2 
in our 54" beam, giving an integrated line intensity 
of 60 K-km s _1 . Given SPIFI's velocity resolution 
of 61 km s" 1 , the intrinsic linewidth is about 50 
km s" 1 . The linewidth and the velocity centroid 
(-35 km s _1 ) are in good agreement with those of 
the radio recombination lines (Brooks, Storey, & 
Whiteoak 2001), but the line is somewhat broader 
and blueward of the neutral gas lines such as CI 
( 3 Pi -^ 3 P ) and CO(4^3) (Zhang et al. 2001), 
consistent with an origin for the [Nil] line in the 
ionized medium. 

4. Discussion 

The Carina Nebula is well-studied morphologi- 
cally and spectroscopically from the X-ray to radio 
regime. Most relevant to the present work is the 
far-infrared spectroscopic mapping obtained with 
the ISO satellite. Mizutani, Onaka and Shibai 
(2002) obtained complete (43 to 197 /im) LWS 
low-resolving-power (R ~ 200) spectra of the Ca- 
rina Nebula. A rectangular region was mapped 
from 1 = 287.°0 to 287.°65 and b = -0.°78 to -0.°23, 
which covers the prominent Carina I and Carina II 
HII regions. The LWS beam size was ~ 70 to 80", 
but the raster scans were coarsely sampled so that 
data were obtained on a square grid with a spac- 
ing of 3'. The prominent far-IR fine-structure lines 
of [OIII], [OI], [CII], [NIII], and [Nil] 122 /an arc 
readily detected at nearly every position mapped. 
The authors use the line ratios to build a model 
for the ionized gas regions in the nebula. In par- 
ticular, the ratio of the 52 to 88 ^im [OIII] lines is 
sensitive to electron density for densities between 
100 and 3000 cm' 3 (cf. Melnick et al. 1979). They 
find two distinct components to the electron den- 
sity: a high-density (n e ~ 100 — 350 cm" 3 ) compo- 
nent at the Carina I and Carina II HII regions, and 
an extended low-density (n e < 100 cm" 3 ) compo- 
nent detectable over the entire ~ 30 pc diameter 
region mapped. Mizutani et al. (2004) find a cor- 
relation between the observed [CII] and [OI] 63 
fim lines from which they estimate that 80% of 



the [CII] line emission comes from PDRs in the 
Carina Nebula as opposed to low-density ionized 
gas regions. 

4.1. The [Nil] Line Ratio 

The [Nil] lines provide a distinct probe of the 
ionization structure and gas density in HII re- 
gions. Formation of ++ requires energetic 35.1 
eV photons, so that the [OIII] lines trace HII re- 
gions formed by very early-type (07) stars. Lower 
energy photons (14.53 eV) can form N + , so that 
these lines arise from HII regions formed by the 
softer UV radiation of late O-type or early B-type 
stars, or from the lower ionization "outskirts" of 
HII regions formed by earlier-type stars. Since 
++ and N ++ have similar ionization potentials 
(35.1 and 29.6 eV, respectively), an ++ zone 
contains mainly N ++ rather than N + . The twin 
[Nil] lines have significantly lower critical densities 
(n c „ t ~ 44 cm" 3 and 293 cm" 3 (8000 K) for the 
205 and 122 /im lines, respectively) for thermaliza- 
tion than the [OIII] lines (n cr i t ~ 510 cm -3 and 
3600 cm" 3 (8000 K) for the 88 and 52 fim lines, 
respectively), so that they make better probes of 
lower-density gas. 

Unfortunately, the undersamplcd ISO/LWS 
mapping of Mizutani et al. (2002) does not di- 
rectly sample the position where we obtained our 
spectrum displayed in Figure 1. The nearest po- 
sition is at 1 = 287.°550, b = -0.°583 (RA = 
10 h 44 m 54 s .l, Dec = -59°37'17" (J2000)), or about 
1/27 (more than a full beam) offset from our peak 
line flux position. However, this position is cov- 
ered by other pixels in the array, for which we 
obtain an integrated line flux of 21.4 K- km s" 1 , 
or 3.6 ± 0.7 x 10" 15 W m" 2 in our 54" beam. 
Using the ISO data archive , we obtain a 122 /im 
[Nil] line flux of 1.66 xlO" 14 W m" 2 , or mak- 
ing the extended source correction, a flux of 1.12 
xlO" 14 W m" 2 in the 78.2" LW2 ISO/LWS detec- 
tor beam (Gry et al. 2003). Assuming a uniform 
intensity source, and correcting for the difference 
in beam size, the ISO flux is 5.33 x 10" 15 W m" 2 
referred to our 54" beam, so that the 122 (im line 
is 1.5 times brighter than the 205 /im line. As- 
suming 30% calibration uncertainty for each line, 
the ratio is 1.5 ± 0.64. 

We have calculated the [Nil] 122 /zm to 205 /im 
line intensity ratio as a function of gas density as- 
suming electron impact excitation and using the 



4 



collision strengths from Hudson and Bell (2004) 
scaled to an assumed electron temperature of 8000 
K (Figure 2). The observed ratio of 1.5 indicates a 
very modest gas density: ~ 32 cm" 3 . Even allow- 
ing for a 30% calibration uncertainty, it is clear 
that the [Nil] lines (tracing low ionization gas) 
arise from a low-density medium: 20 < n e < 60 
cm" 3 . Therefore, the "halo" of low density ionized 
gas discovered in its [OIII] line emission (Mizutani 
et al. 2002) also contains gas in lower ionization 
states. This low ionization, diffuse gas is similar 
to the warm ionized medium in the Galaxy as a 
whole. 

4.2. The Fraction of [CII] from Ionized 
Gas 

Since the critical densities for electron impact 
excitation of the 158 ^m [CII] line and the 205 /im 
[Nil] line are very similar, to a good approxima- 
tion the line ratio in ionized gas regions is depen- 
dent only on the relative abundance of C + and N + 
within the HII region. Using the collision rates 
for exciting the ground-state levels of C + from 
Blum and Pradhan (1992), we have calculated the 
expected ratio of the two lines as a function of 
density and present it in Figure 2. The temper- 
ature dependence is quite small, as the levels are 
only 91 and 70 K above ground, respectively — 
small compared with the temperature (8000 K) of 
an HII region. For the calculation, we take the 
gas-phase abundances of n(C + )/n e = 1.4 x 10" 4 
and n(N + )/n e = 7.9 x 10~ 5 (Savage and Sembach 
2004) and correct for the fraction of the two el- 
ements expected in the first ionization state us- 
ing the HII region models of Rubin (1985). The 
[CII] /[Nil] line ratio from ionized gas is 3.1 at low 
densities, and 4.3 at high densities, and dips to be- 
low 3 between densities of 20 and 100 cm" 3 since 
the upper J=2 level of the N+ begins to be signif- 
icantly populated at densities above 20 cm" 3 . 

From the ISO archive, we obtain a [CII] flux 
of 1.01 x 10" 13 W m" 2 . Making the extended 
source correction, this corresponds to 5.5 x 10" 14 
W m" 2 within the 70" (Gry et al. 2003) LW4 
beam. Again, assuming a uniform source the flux 
within our 54" beam would be 3.3 x 10" 14 Win" 2 . 
Allowing for a 30% calibration uncertainty in each 
line, the [CII] /[Nil] 205 line ratio is - 9.2 ± 
3.9 to 1. Inspecting Figure 2 for a gas density of 32 
cm" 3 , if both lines arise from the ionized gas, we 



would expect the line ratio to be 2.44. Therefore, 
27% (with an upper bound of 46% and a lower 
bound of 19%) of the observed [CII] flux arises 
from the ionized medium, and the remaining 73% 
must have its origin in the neutral ISM. These re- 
sults support and underpin prior work which con- 
tends that most of the observed [CII] line emission 
from Galactic star- forming regions, the Galaxy as 
a whole, and from external galaxies arises in warm 
dense photodissociation regions on the surfaces of 
molecular gas clouds. 
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Fig. 1. — [Nil] 205 (im spectrum obtained using 
SPIFI on AST/RO. These data were obtained at 
the peak of the radio continuum contours of the 
Carina II HII region (Retalack, 1983). Velocity 
resolution is 60 km s _ l. The dashed line shows 
the intrinsic instrumental profile. The data were 
smoothed with a Hann window. 
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Fig. 2.— [Nil] 122 nm to 205 yum line intensity 
ratio as a function of ionized gas density, n e . (solid 
line). The open circle marks our observed line 
ratio of 1.5 ±0.64. The ratio of the [CII] 158,umto 
[Nil] 205 /im lines as a function of n e (dashed line). 
The filled circle marks the observed line ratio of 
9.2 ± 3.9. 
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